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Semi-metals are characterized by nodal band structures that give rise to exotic electronic proper¬ 
ties. The stability of Dirac semi-metals, such as graphene in two spatial dimensions (2D), requires 
the presence of lattice symmetries, while akin to the surface states of topological band insulators, 
Weyl semi-metals in three spatial dimensions (3D) are protected by band topology. Here we show 
that in the bulk of topological band insulators, self-organized topologically protected semi-metals 
can emerge along a grain boundary, a ubiquitous extended lattice defect in any crystalline material. 
In addition to experimentally accessible electronic transport measurements, these states exhibit 
valley anomaly in 2D influencing edge spin transport, whereas in 3D they appear as graphene-like 
states that may exhibit an odd-integer quantum Hall effect. The general mechanism underlying 
these novel semi-metals - the hybridization of spinon modes bound to the grain boundary - sug¬ 
gests that topological semi-metals can emerge in any topological material where lattice dislocations 
bind localized topological modes. 


Graphene [T] and topological band insulators (TBIs) 
II [3] show exotic electronic transport properties that 
have motivated the search for other materials exhibit¬ 
ing similar semi-metallic features. Semi-metals are de¬ 
scribed by electronic band structures where the bands 
touch at isolated points or lines in the Brillouin zone 
(BZ). In graphene, a 2D honeycomb lattice of carbon 
atoms, or in Dirac semi-metals in 3D, the bulk hosts a 
pair of pseudorelativistic gapless Dirac fermions, while 
the surface states of TBIs feature in general gapless Weyl 
fermions - chiral massless particles extensively studied in 
high-energy physics for the description of neutrinos. Re¬ 
cently, the latter have been discovered also in the bulk 
of 3D materials known as Weyl semi-metals HHl]. While 
the energy spectra in all cases resemble each other, the 
stability of their band structures has a dramatically dif¬ 
ferent origin. In Dirac semi-metals the stability of the 
Fermi surface relies on lattice symmetries, while Weyl 
semi-metals and surface states of TBIs are protected by 
the topology of the bulk band structure. Therefore, it is 
of both fundamental and practical importance to answer 
the following question: Can topologically protected Weyl 
fermions also appear in the bulk of lower dimensional 
systems? 

We here provide an affirmative answer to this question 
by showing that grain boundaries (GB) - ubiquitous crys¬ 
tal defects in real materials that are usually considered as 
detrimental for their properties - can host time-reversal 
symmetry (TRS) protected topological semi-metals. GBs 
arise at the interface of two crystal regions (grains) whose 
lattice vectors are misaligned by an angle 0, as illustrated 
in Fig. lA. For small opening angles a GB can be viewed 
as lattice dislocations described by Burgers vector b ar¬ 
ranged on an array of spacing d = |b|/(tan 0). While GBs 


are usually considered as unwanted disorder, they have 
also been used experimentally as probes of the supercon¬ 
ducting state in high-temperature superconductors it 
HU]. Recently, they have also been suggested for engi¬ 
neering thermoelectric devices m and for experimen¬ 
tally tuning the surface states in a 3D TBI [12] . Our main 
result is that extended GB lattice defects can host stable 
self-organized states of matter. We show that GBs in the 
bulk of 2D and 3D TBIs can realize stable TRS protected 
ID and 2D semi-metals, respectively, which, in contrast 
to Dirac fermions in graphene, do not exhibit pseudospin 
degeneracy. These ’’halved” graphene-like states can be 
experimentally observed by measuring their characteris¬ 
tic conductance through otherwise insulating bulk. They 
are also intimately connected to the TBI surface states 
and can influence surface transport: ID GBs in the bulk 
of a 2D TBIs exhibit a valley anomaly that under an ap¬ 
plied electric field results in a helical imbalance of the 
TBI edge states on the two grains connected by the GB. 
Furthermore, 2D GBs in the bulk of 3D TBIs may exhibit 
an odd-integer quantum Hall effect. 


I. SEMI-METALLIC STATES ON GRAIN 
BOUNDARY LATTICE DEFECTS 

The physics of GBs can be derived from their elemen¬ 
tary building blocks - the lattice dislocations. It has been 
shown that when a 2D TBI is in a so called translationally 
active topological phase [T3| , where the band topology is 
characterized by an odd number of band inversions at 
TRS momenta other than the F-point in the BZ, a lat¬ 
tice dislocation acts as an effective magnetic tt-Aux and 
binds a single Kramers pair of zero energy spinon modes 
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FIG. 1. Graphene-like semi-metal on a ID grain boundary in a translationally active 2D TBI (the M-phase of the BHZ model 
on a square lattice, see Appendix [A| . (A) Schematic illustration of a GB. The coordination discrepancy due to the angular 
mismatch 0 of lattice basis vectors results in an effective array of dislocations, each described by Burgers vector b and marked 
by a ’T’ symbol, of spacing d = |b|/ tan^. (B) Real space numerical tight-binding calculation of the low-energy states in the 
presence of a GB. The spinons bound to the dislocations hybridize with tunneling amplitude t, which gives rise to an extended 
ID state along the GB. The radii of the circles indicate the amplitude of a wave function, associated with the node at kx — 'k jd 
of panel C, while the colors indicate the phase. (G) The characteristic mid-gap bow-tie dispersion of the hybridized spinon 
bands (orange) and bulk bands (grey), corresponding to GB opening angle 0 = 18.4° for which isolated dislocations along the 
GB are well dehned, as the function of the momentum 0 < < 27r/d along the GB. (D) The same plot for the maximal 

opening angle 0 = 45°. While the picture of isolated dislocations breaks down and the simple nearest-neighbour hybridization 
based sinusoidal dispersion is lost, the defining TRS semi-metal dispersion with the nodes at the TRS momenta survives. The 
data are for a 30 x 60 site or larger system with periodic boundary conditions. 


localized at the core [HHni. These modes are 2D ana¬ 
logues of the famous spin-charge separated Jackiw-Rebbi 
soliton states [iHlIlS] realized in the Su-Schrieffer-Heeger 
(SSH) model of polyacetylene [20l[2T]. The same mech¬ 
anism applies also to 3D TBIs in translationally active 
phases, which generalize the characterization by weak in¬ 
dices [Tang. The natural lattice defect is a ID disloca¬ 
tion line that behaves as tt-Aux tube and binds counter- 
propagating helical modes akin to edge states in 2D TBIs 
[23l[24]. A complete catalogue of modes bound to dis¬ 
locations in generic 3D TBIs is given by the K — b — t- 
rule, that relates the spectrum of the surface states to 
the number of spinon modes bound to a dislocation line 
piercing the surface m- 

When the spinon zero modes form an array, such as 
along a GB, one expects them to hybridize into an ex¬ 
tended state. In the basis of two spinon modes, this 
state should be described by a two-band model h{q) = 


ho{q) + h(g) a. As a act in the spinon basis, TRS re¬ 
quires hQ{—q) = ho{q) and h(—g) = —h{q). Thus based 
on these very general considerations only, one expects a 
semi-metal with two TRS protected nodes at the TRS 
momenta q = 0 and q = tt. These nodes are degener¬ 
ate when ho{q) either vanishes or is a constant, which, 
as we show below, is determined by the presence of bulk 
inversion symmetry. 

To verify this prediction, we have carried out numer¬ 
ics for GB defects in the translationally active phases of 
BHZ tight-binding models for topological insulators with 
both TRS and inversion symmetry [25] (see Appendix [A| 
for details). Indeed, Figuresandshow that GBs sup¬ 
port propagating states along their core. In the energy 
spectrum these states appear as emergent mid-gap ’bow- 
tie’ bands that, reminiscent of the valleys in graphene, 
exhibit the expected two degenerate nodes at the two dis¬ 
tinct TRS momenta. We have verified that these bands 









































appear throughout the translationally active phases and 
persist for all GB angles including the maximal opening 
angle of ^ = 45° close to which a GB can no longer be 
approximated by an array of isolated dislocations. Thus 
when the chemical potential of the parent TBI is at the 
nodes, a self-organized graphene-like semi-metal emerges. 


II. EFFECTIVE MODELS FOR SEMI-METALS 
ON GRAIN BOUNDARIES 

Having numerically established the existence of ex¬ 
tended states on grain boundaries, we now turn to con¬ 
struct effective models for these semi-metals. The micro¬ 
scopic mechanism underlying their emergence is the hy¬ 
bridization of topological spinon modes bound to disloca¬ 
tions. We derive the hybridization induced dynamics by 
employing the construction for dislocation bound spinons 
that relies only on the existence of edge states originat¬ 
ing from the band inversion momenta and their sym¬ 
metries [23]. The edge states relevant to construct the 
spinon modes can be identified using the K — h — t-rule 

E31- 

Gonsider a system of two coupled helical edges along 
^-axis described by the Hamiltonian 

Ho = vkyCF^ iJi^ + ( 1 ) 

where acts in the spin space, while = ±1 denotes 
the two edges. In general, tunnelling of magnitude m 
gaps out the edge states at the interface of the two TBIs. 
However, when an additional row of atoms for ^ > 0 is 
inserted between the edges, or equivalently a dislocation 
with Burgers vector b = is created at ^ = 0, the 
mass term for ^ > 0 becomes me* |23l|2lj. In 

a translationally active phase ^ • b = tt (mod 27r) 

and thus the sign of the edge state tunneling becomes 
^-dependent. The system is then equivalent to the ID 
continuum SSH model with a mass domain wall EHIET]. 
In a TRS system the domain wall binds a Kramers pair 
of localized spinons, which for the Hamiltonian Hq live in 
a subspace described by the projector P = (1 -hcr^/i^)/2. 

To model a GB in a 2D translationally active topolog¬ 
ical material, we consider a dislocation array of spacing 
d along x-axis. Employing the construction described 
above, the adjacent edge states then couple by a transla¬ 
tionally and time-reversal invariant Hamiltonian 

Hgb = t{cos kxii^ + sin k^lJ^), (2) 

where kx spans the reduced BZ 0 > kx > 27r/d. The 
coupling strength t follows from the overlap of the edge 
states, that is proportional to the overlap of the localized 
spinon wavefunctions (see Appendix]^, and thus scales 
as t ^ implying that it is stronger for large opening 
angle GBs. Projecting this coupling into the spinon sub¬ 
space, consistent with our prediction, we find that the 
spinons acquire dynamics described by 


( 3 ) 


where is an effective spin operator in the spinon sub¬ 
space (see Appendix]^. Gomparing this form to our nu¬ 
merics, Figure shows that this simple expression in¬ 
deed captures the defining nodal structure of the emer¬ 
gent mid-gap band for low opening angle GBs in the 
presence of spin-rotational and inversion symmetries (the 
large opening angle GBs require couplings beyond near¬ 
est neighbour). 

In Appendix we have considered in detail the gen¬ 
eral case of breaking all the symmetries of the parent 
BHZ model. In the effective model Sz conservation 
breaking via a Rashba term is modelled with Hr = 
akyCF^^ inversion breaking with Hj = (m^ 

can be either constant or proportional to cos kx) and TRS 
breaking via the Zeeman terms Hr = In the 

presence of the Rashba term the projector to the spinon 
subspace is given by 
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and the effective Hamiltonian H 2 D = Pa{HGB + Hj + 
Hr)P a becomes 


H2D 


rriya + rrizV 


vhz + ahy 


t sin k 


CF^. (5) 


This expression shows that the key property determin¬ 
ing the response to TRS and inversion breaking is the 
spin texture of the edge states: Only perturbations that 
couple to a spin orientation present in the edge states ap¬ 
pear in the effective model (e.g., if Sz is conserved, only 
perturbations proportional to appear in the effective 
Hamiltonian). 

This effective picture is fully consistent with our nu¬ 
merics on the stability of the GB state (see Appendix 
1^ for the details of the numerical calculations). First, 
Sz conservation breaking Rashba coupling that is not 
strong enough to close the bulk band gap preserves the 
semi-metallic nodes. Second, terms breaking the bulk in¬ 
version symmetry in general shift the nodes to different 
energies, but, in contrast to interfaces of TBIs with dif¬ 
ferent velocities [26], can not gap them out. Third, TRS 
breaking terms gap the spectrum, but only if they do not 
anti-commute with Po, and only once their magnitude 
is comparable to the bandwidth 2t resulting from the 
hybridization. Furthermore, we have numerically found 
that no moderate random disorder can open up a gap 
at the nodes since the topological edge states underly¬ 
ing the hybridizing dislocation modes persist as long as 
the disorder does not drive the bulk out of the trans¬ 
lationally active phase. This also implies that broken 
translational invariance along the GB, that can arise due 
to lattice reconstruction leading to disordered dislocation 
positions, does not pose a fundamental obstacle for the 
stability of the nodal band structure. Similar to chemical 
potential disorder that deforms dislocation mode wave- 
functions, random bends along the GB result in random 
tunneling couplings t. In terms of the low-energy theory 


P{HGB)P = tsmkxdy 
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FIG. 2. Graphene-like semi-metal localized on a 2D grain boundary in a translationally active 3D topological band insulator 
(the i?-phase of the BHZ model on a cubic lattice, see Appendix]^ for details). (A) Schematic illustration of the 2D GB, which 
now realizes a sheet of parallel ID edge dislocation lines that extend along the grain boundary. Each edge dislocation hosts a 
pair of propagating helical modes localized along the dislocation cores that cross at kz = tt. (B) On a GB of opening angle 18.4°, 
the hybridization of the helical modes results in the semi-metal band structure with two anisotropic pseudo-relativistic fermions 
at (/Ccc, kz) — (0,7r) and (vr/d, tt). Along the GB we recover the same mid-gap bow-tie dispersion as in the 2D case (shown top 
right for kz = tt), while along the edge dislocations the hybridized modes (orange) still flow into the bulk bands (gray) as a 
function of kz (bottom right for kx — 7r/d). The data are for a 60 x 60 x 90 system with periodic boundary conditions. 


around the nodes, this gives rise to a random gauge field 
that only shifts the cones [27j. As long as this shift is 
smaller than the separation 7r/d of the cones in the BZ, 
the semi-metallic behavior is stable. We have numeri¬ 
cally verified this argument by considering random dis¬ 
location positions along the GB and found qualitatively 
similar stability to the case of chemical potential disor¬ 
der. The self-organized semi-metal on the GB thus shares 
the topological stability of the edge states of the parent 
state: The nodes are degenerate in energy only if the 
bulk inversion symmetry is intact, but TRS is sufficient 
to protect the nodes themselves (no additional spatial 
symmetries with respect to the GB need to be assumed). 

This mechanism generalizes straightforwardly to 3D 
translationally active TBIs, where a GB consists of a 2D 
sheet of parallel ID dislocation lines, as illustrated in 
Fig. |2}\. In a translationally active phase each disloca¬ 
tion binds a Kramers pair of helical modes [23l[24]. To 
derive a minimal effective model for their hybridization, 
we employ again the surface state construction with a 
coupling between adjacent surfaces. In the geometry of 
this is described by the Hamiltonian 

H = v{ky(jy + + m{y)ix^ + Hgb- (6) 

Similar to the 2D case, we project this Hamiltonian into 
the spinon subspace near kz = 0, which in the presence of 
TRS breaking Zeeman terms Hb and inversion breaking 


terms Hj gives the general minimal model (see Appendix 

H 

Hsd = rriy ^ vkzd^ + {hy + t sin k^) (7) 

This effective theory describes two linearly dispersing 
cones at the two TRS momenta {k^^kz) = (0,7r) and 
(yr/d, tt) with anisotropic velocities v along and t per¬ 
pendicular to the dislocation lines, respectively, consis¬ 
tent with Fig. [^. This anisotropy is reduced for larger 
opening angles 0 due to an enhanced overlap between 
the helical modes whose hybridization underlies the self- 
organized GB state. The two cones separated in the BZ 
are degenerate in energy if the bulk inversion symme¬ 
try with respect to the GB plane is preserved {rriy = 0). 
On the other hand, if inversion is broken, the cones ap¬ 
pear at different energies. As with the 2D case, we have 
numerically verified that the semi-metal is stable in the 
presence of moderate disorder and that the cones can 
only be gapped out with a TRS breaking Zeeman terms 
normal to the GB plane (see Appendix [D|). 

III. EXPERIMENTAL CONSEQUENCES 

Having established the topological stability of the 
emergent semi-metals on GBs, we now turn to their ex¬ 
perimental signatures. To detect the semi-metallic state 
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FIG. 3. The experimental setups for identifying the signatures of the graphene-like semi-metals. (A) The one-dimensional 
bow-tie dispersion of the spinon semi-metal on a ID GB implies a parity anomaly per spin component when an electric field 
E is applied along the grain boundary. The arrows indicate the shift in the spectrum from one valley to the other: When Sz 
is conserved, valley vi accumulates an excess of spin down (red), while valley V 2 accumulates excess spin up (blue). These 
valleys can be associated with two co-existing channels that connect the helical edge states from the opposing surfaces. When 
E is applied along the GB, a current for both spin orientations is driven parallel to it. At the GB termination points this 
current is predicted to flow into the edge states that propagate to opposite directions resulting in a doubled spin Hall effect-like 
helical imbalance of the edge states on the two grains. Measuring the current imbalance I is the hallmark signature of the 
valley anomaly exhibited by the spinon semi-metal. (B) The spinon semi-metal on a 2D GB may feature an odd-integer Hall 
effect with transverse conductivity Gxz — (2n + l)e^/h in the presence of the perpendicular magnetic field B. In the absence 
of external fields, another signature is provided by the diagonal ballistic optical conductivity of azz = {ttI 4:)e^ which is half 
that of graphene. 


on the ID GB inside a 2D TBI, one can carry out a di¬ 
rect two terminal transport measurement analogous to 
the one used to detect edge states in a 2D quantum spin 
Hall insulator [28]. When an electric field E is applied 
along the GB with leads attached to the GB ends, due 
to the two cones one should observe conductance of 

a = jh, ( 8 ) 

i.e. twice the value measured for the QSH edge states. A 
more dramatic consequence of the helical bow-tie disper¬ 
sion along a ID GB is the existence of the parity anomaly 
m for each helical band. In the presence of two chiral 
cones, it gives rise to a valley anomaly that influences 
the edge transport. As illustrated in Fig. [^, an electric 
field E along the GB generates excess helicity of opposite 
orientations at the two valleys that results in a net cur¬ 
rent flowing along the applied field. Associating the two 
valleys with two co-existing channels through which the 
helical edge states can flow from one GB termination sur¬ 
face to the other, the valley anomaly results in a helical 


imbalance of the edge currents at the two grains on each 
side of the GB. This imbalance is proportional to the hy¬ 
bridization strength t and it represents a hallmark trans¬ 
port signature of valley anomaly. To detect it, one can 
carry out a two-terminal edge transport measurement il¬ 
lustrated in Fig. [^. When the two edges connected by 
the GB are biased by voltage a net current I ^ 
is expected as a consequence of unequal helical currents. 
Recent experimental observation of semi-metallic trans¬ 
port on domain walls [30l|3T], such as in bi-layer graphene 
j32H35] , is highly encouraging that an experiment of this 
kind can also be carried out for GBs. In materials where 
spin is a good quantum number, i.e. spin-orbit coupling 
is negligible, the helical imbalance due to applied electric 
field translates into a spin imbalance of the edge currents 
on the two grains. In such cases the valley anomaly could 
also be detected by measuring the magnetic moment due 
to spin imbalance |36| or using Kerr rotation microscopy 
|37l[38]. 

While the cones of the graphene-like semi-metal on 







































a 2D GB can also be viewed as two co-existing chan¬ 
nels through which the surface states of the 3D TBI 
can propagate between the surfaces connected by the 
GB, this state is experimentally most conveniently de¬ 
tected via two distinct ” half-graphene”-like transport sig¬ 
natures. First, the measured ballistic optical conductiv¬ 
ity is Gzz — I which is half the value measured 
in graphene [39]. Second, when a magnetic field is ap¬ 
plied perpendicular to the GB in the setup shown in Fig. 
[^, the non-degeneracy of each cone implies a contribu¬ 
tion of {n ll2)e^ jh to the Hall conductivity, which in 
turn implies odd-integer quantum Hall effect with total 
Hall conductivity 

a^z = (2n + l)e^/h, (9) 

when the cones are degenerate in energy (bulk inversion 
symmetry is present). As there is always a known contri¬ 
bution from the surface states, it is in principle possible to 
extract the signal associated with the Hall conductance 
arising solely from the GB semi-metal. The emergent chi¬ 
ral symmetry of the effective model 0 also suggests the 
existence of edge states on the ID GB edges on the 2D 
surfaces (see Appendix H- These can contribute addi¬ 
tional Fermi arc-like features to the surface states of the 
parent 3D TBI, which may be detectable via zero fre¬ 
quency optical conductivity measurements m , or give 
rise to GB edge transport when the cones occur at dif¬ 
ferent energies due to broken inversion symmetry m- 


IV. DISCUSSION AND OUTLOOK 

We have shown that TRS protected semi-metals - he¬ 
lical wires exhibiting a valley anomaly in 2D and an 
anisotropic graphene-like state showing an odd-integer 
quantum Hall effect in 3D - can emerge on grain bound¬ 
aries in TBIs. These states are self-organized and enjoy 
the same topological stability as the TBI edge states. 
The only requirement for their emergence is the TBI to 
be of a translationally active type where lattice disloca¬ 
tions bind a Kramers pair of localized modes. As grain 
boundaries occur naturally in crystalline samples, or they 
can be experimentally manufactured [42] , the challenge is 
to identify translationally active materials. A convenient 
guide to candidates is given by the space group classi¬ 
fication of TBIs [13]. For modelling purposes we used 
the M- and R-phases of the 2D square and 3D cubic 
lattice, respectively, with the latter having a potential 
realization in electron-doped BaBi02 [43]. Nonetheless, 
our mechanism is completely general and the outlined 
GB physics should therefore also appear in the already 
experimentally verified translationally active materials. 
This rapidly expanding list includes Bia,Sbi_a^ with band 
inversions at the L-points [44], the topological Kondo 
insulator SmBe with band inversions at three X points 
[45] [46] and the recently discovered bismuth iodide com¬ 
pounds that feature a transitionally active phase with 


band inversions at the Y points 07 ]. Topological semi¬ 
metals are expected also on GBs in the topological crys¬ 
talline Sn-based compounds [48ll5Q] . that feature band 
inversions at the L points as long as the protecting sym¬ 
metry is respected iniiiii, as well as in Bismuth bilayers 
with a band inversion at the M point, where a GB state 
has been identified in a recent ab initio study [51]. 

Grain boundaries provide a natural setting for arrays 
of localized topological solitons to hybridize in any crys¬ 
talline materials. The emergent extended states depend 
on the topology and the symmetry class of the parent 
state. Our general mechanism is readily applicable to 
model the distinct GB states that may be realized in 
topological states of matter with symmetries different 
than considered here. For instance, lattice dislocations in 
Sr 2 Ru 04 , a candidate for chiral p-wave superconductor, 
can host Majorana modes [52], which are expected to hy¬ 
bridize into novel superconducting states [53] . A full clas¬ 
sification of extended states on GBs in crystalline topo¬ 
logical states of matter provides for a fascinating subject 
of future work. Finally, our results demonstrate the po¬ 
tential of extended lattice defects in exploring manifes¬ 
tations of electronic band topology beyond by now well 
understood surface states. 
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Appendix A: Tight-binding BHZ models for 
numerics 

1. The 2D BHZ model 


The 2D BHZ model is defined on a square lattice m 
with two spin degenerate \s) and \px + ipy) type orbitals 
on every latttice site. In natural units h = c = e = the 
model is defined by the nearest neighbour tight-binding 
Hamiltonian 


Utb = + H.C.) + ^ (Al) 

r,S r 

Here denote the vectors con¬ 
necting the nearest neighbor sites and = 

( 5 ^(r),p>^(r), s^(r),p^(r)) annihilates the s and p type or¬ 
bitals at site r. The tunneling of the spin up and spin 
down electrons is given by 


r i* _ rji 

5,u ~ 


f As ts/2\ 

\t'^/2 Ap J 


that describes inter-orbital tunneling ts = —iex.p{i(ps) 
and = —iexp{—iips) whose phase is given by the po¬ 
lar angle cps of the vector and intra-orbital tunnel¬ 
ing of magnitude Ag/p = AB. The on-site energies are 
parametrized as /i = (M — 4:B)tz ( 8 ) ctq, where the Pauli 
matrices r and a act in the orbital and spin space, re¬ 
spectively. Here, tq and do are the 2 x 2 identity matri¬ 
ces. By performing a Fourier transform, the Hamiltonian 
(Al) assumes the block-diagonal form 


Vt fH{^) 0 


Utb = 


(A 2 ) 


as being translationally active. In this phase the edge 
states appear always at the projection of the M-point to 
the surface BZ, i.e. they cross at /c = tt on the edge BZ. 

The 2D BHZ model has additional symmetries besides 
TRS. Explicitly, all the symmetries comprise: 

• Time-reversal symmetry TH(k)T~^ = H{—k) rep¬ 
resented by T = i{(Ty (8) To)A satisfying = — 1, 
where and K denotes complex conjugation. 

• Particle-hole symmetry PH(k)P~^ = —H{—k) 
represented by P = (do 0Tx)K satisfying P^ = 1 . 

• Chiral symmetry CH{k)C~^ = —P[{k) represented 
by C = PT = i(jy ( 8 ) Tx . 

• Spin-rotation symmetry SzH{k)S~^ = H{k) rep¬ 
resented by Sz = C) To. 

• C '4 crystalline point group symmetry in the x-y- 
plane represented by P = 

that maps kx^y Akx^y^ t^’^ ^ and 

It gives rise to an inversion sym¬ 
metry IH{k)I~^ = P[{—k) satisfying P = —1 rep¬ 
resented by / = = ia^^Tz, which can be further 

decomposed into inversions along x- and y-axes as 
/ = Ixly, where Ix = icr^P and ly = obey 

IiH{ki)Ip = H{-ki). 

The Sz symmetry implies that the model can be viewed 
as two copies of opposite chirality Chern insulators, one 
for each spin component. This conservation is non¬ 
generic and is in general broken by spin-orbit coupling. 
For the BHZ model specific to HgTe quantum wells, the 
orbital dependent Rashba form reads 


The Hamiltonians for each spin component can be de¬ 
composed as P(k) = T • d(k), where the vector d(k) 
has the components dx^yi)^) = Psm{kx,y) and dz = 
M — 2B{2 — cos{kx) — cos{ky)). The spin up and down 
blocks are related by the time-reversal symmetry repre¬ 
sented by an antiunitary operator T = tq (8) with 

K as the complex conjugation. The energy spectrum for 
each spin component is given by P(k) = ±|d|, while the 
spectrum of the full Hamiltonian is doubly degenerate. 

The dispersion P(k) is gapped except for the values 
M/B = 0,4 or 8 , where the gap closes at the F (0,0), 
X ( 7 r, 0 ) and Y ( 0 , 7 r), or M ( 7 r, 7 r) points of the Bril- 
louin zone, respectively. When 0 < M/B < 4 the sys¬ 
tem is in a topological F-phase, while for 4 < M/B < 8 
the system is in a topological M-phase. For other val¬ 
ues of M/B the system is topologically trivial and does 
not have helical edge states. The F- and M-phases are 
characterized by the same Z 2 invariant and both exhibit 
helical edge states. They are distinguished by lattice de¬ 
fects that break translational symmetry: In the F-phase 
a lattice dislocation does not bind localized modes, while 
in the M-phase they act as 7 r-fluxes and bind localized 
zero energy modes m- Hence, we refer to the M-phase 


Ts ^ Ts P + 'T~z)^z • (cr X ^). {AS) 

This term breaks also the accidental P and C symme¬ 
tries, leaving the model with only TRS and the C 4 point 
group symmetry. The latter can be further broken, while 
preserving TRS, by terms of the form Hi = 
with non-zero rrii breaking the inversions li. The coef¬ 
ficients rrii are constant if they originate from local po¬ 
tentials, or proportional to cos ki when they result from 
tunneling processes. 


2. The 3D BHZ model 

The 2D BHZ model on a square lattice can be gener¬ 
alized to a 3D cubic lattice. The key difference is that in 
3D the model can no longer be viewed as two decoupled 
Chern insulators of opposite chirality, and an Sz con¬ 
servation breaking spin-orbit coupling is needed for the 
model to exhibit topological phases. Working directly 
in the momentum space, the tight-binding model con¬ 
tains two spin degenerate orbitals, which are in general 



described by a Hamiltonian of the form 

H = e(k)l + da{^)ja + dal3l(xl3i 

oc a/3 

where 7 q, are the five Dirac matrices obeying the Clifford 
algebra { 7 a, 7/3} = 7 a /3 are the ten commutators 

7 a /3 = ^[ 7 a, 7 / 3 ]- Specifically, we take the following rep¬ 
resentation for the 7 -matrices 


7o — cTo < 


73 = 


7l= l2=cry 


(A4) 


75 = -70717273 = cro 0 Ty , 


with a and r being again the standard Pauli matrices 
acting in the spin and orbital space, respectively. Tak¬ 
ing into account only nearest-neighbor terms and setting 
the lattice constant to identity, the BHZ model in 3D 
is described by the Hamiltonian (A2) where d^is = 0, 
do = M — 2B{3 — coskx — cosky — coskz) and di^ 2,3 = 
sin kx^y^z- As with the 2D model, B is again the magni¬ 
tude of the intra-orbital hopping and M is the difference 
of the onsite energies between the two orbitals. How¬ 
ever, one should keep in mind that while this Hamiltonian 
bears formal similarity to the 2 D one (A 2 ), the explicit 
spin-orbit coupling is included in the 7 -matrices that de¬ 
scribe spin-flipping inter-orbital tunneling, whose magni¬ 
tude has been set to unity. The phase diagram exhibits 
now topological phases for 0 < M/B < 4, 4 < M/B < 8 
and 8 < M/B < 12. Only the last, the i?-phase with 
band inversion only at the i?-point (7r,7r,7r), is a trans- 
lationally active phase where the surface states always 
cross at (tt, tt) on the surface BZs m- 

This model has also additional symmetries originating 
from the orbital space. In addition to TRS, the model 
has the global particle-hole and chiral symmetries de¬ 
scribed by P = Ka^ry and C = ir^, respectively. The 
cubic lattice point group symmetry gives rise to inver¬ 
sion I = which can be again decomposed into inver¬ 
sion li = normal to each of the three axes. The 

P and C can again be broken by Hj = mia^y with¬ 
out breaking TRS, with each non-zero rui breaking again 
also the inversion li. Even in the presence of such terms, 
the 3D model still has a residual particle-hole symmetry 
r^P(k)r^ = — P(—k), which can be further broken by 
moT^ that also breaks all inversion symmetries. 


3. Grain boundaries in tight-binding models 

To simulate a GB in the tight-binding numerics, we 
consider a lattice that consists of two parts with differ¬ 
ent number of rows glued together across the GB defect. 
A lattice with a GB of spacing d with translational in¬ 
variance in ^-direction is created by taking the two halves 
with Lx and Lx/idP 1 ) rows of sites and leaving every 
(d + l)th row unconnected, a representative of which is 
illustrated in Fig. Finally, we corroborated our results 
with an alternative procedure of simply matching two 


perfect crystalline regions under an angle, schematically 
illustrated in Fig. lA of the main text. The resulting 
system with coordinate discrepancy then produced the 
same qualitative results, as illustrated with the 45° re¬ 
sult presented in the main text. 



FIG. 4. Illustration of the grain boundary lattice used in the 
simulations in the particular instance of spacing d = 8. 


Appendix B: Derivation of effective models for the 
GB states 

1. An effective model for the ID GB inside a 2D 
TBI 

We find the effective theory for the GB by coupling 
together adjacent the helical edge states of QSH insu¬ 
lator with a periodic array of trenches with dislocation 
lines inserted, following the construction of Ref [23] for 
solitons bound to dislocations. In the absence of Rashba 
coupling Sz is conserved and two decoupled helical edge 
states along ^-axis are described by the Hamiltonian 
































Hq = vkyCj^ii^ ^ where the a-matrices act in the spin space 
and the = ±1 denotes the two edges. Proximity tun¬ 
neling between the edge states is described by the term 
m/i^, which in general gaps the edge states and merges 
the two QSHE insulators into a single one. When a semi¬ 
infinite row of atoms is inserted for ^ > 0, creating a dis¬ 
location of Burgers vector b = at ^ = 0, the tunneling 
term giving mass to the edge states across the inserted 
row of atoms becomes m [23l [24], where 

Ki are the band inversion momenta from which the edge 
states originate. In other words, in a translationally ac¬ 
tive phase where b • = tt (mod 27r), the mass term 

becomes position dependent such that m{y < 0) > 0, 
while m{y > 0) < 0. The low-energy theory for a single 
trench is then given by 

Hq = vkya^/i^m{y)/i^^ (Bl) 

which is equivalent to a continuum theory for a SSH 
model with mass changing sign at ^ = 0. Such model 
is known to have two bound state solutions at the mass 
domain wall, one for each spin component. The gen¬ 
eral form of these localized solutions is given hy ip = 
ev S with the four-spinor obeying 11 ^( 1 ) = (p. 

This condition defines the subspace where the localized 
soliton solutions live, with the projector onto this sub¬ 
space P = (1 + cr^/i^)/2. Considering a translationally 
invariant array of such bound modes, their hybridization 
can be modelled by coupling the edge states from adja¬ 
cent trenches by introducing 

Hgb = t (cos(g)/i^ -h sin(g)/i^). (B2) 

Projecting this down to the spinon subspace yields the 
minimal effective model PHqbP = tsin(g')d^, where the 
effective spin operator is defined by = jp. 

The effective model Pq + Hqb has additional symme¬ 
tries represented by Sz = P = C = 

and I = ia^y^ symmetries besides TRS. To derive an 
effective model for the fully symmetry broken case, we 
consider edge states with generic spin texture Hq = 
k{v ■ a)y^ m{y)y^, where we have defined the normal¬ 

ized vector V = {vx,Vy^Vz)/\^^\. The spinons live then in 
a subspace defined by the projector P = [1-f (v-cr)/i^]/2, 
which breaks the Sz, P and C symmetries. To con¬ 
sider the effect of breaking the inversion symmetry and 
TRS, we introduce the perturbations y^ 

and Ht = ^i^iCF^, respectively. Here the magnitudes 
hi of the Zeeman terms are constants, while in general 
the inversion-symmetry breaking arises from edge state 
tunneling and thus rrii ^ cos{q). Projecting all these 
down to the spinon subspace gives the effective model 
H 2 d = P{Hgb + P/ + Ht)P, or explicitly 

H 2 D = m • V + (t sing -t- h • v) (B3) 

The spectrum is given by E± = m • v ± (t sin g + h • v), 
which for v = (0, a, v)/\/S v‘^ reduces to the one pre¬ 
sented in the main text. The key observation is that the 


effect of the (spin-rotational symmetry breaking) spin- 
orbit coupling is to change the spin texture of the edge 
states, and that controls which terms perturb the effec¬ 
tive model. Only inversion or TRS breaking perturba¬ 
tions with a component parallel to the spin texture ap¬ 
pear in the effective model, which explains why in our 
numerics only certain perturbations affect the spectrum 
in the full 2D BHZ model. In general, inversion breaking 
terms shift the nodes at the TRS momenta to different 
energies, while TRS breaking terms gap them out once 
their magnitude is of the order of the bandwidth t result¬ 
ing from the hybridization. Other TRS invariant tunnel¬ 
ing terms can also be added to the Hamiltonian, but we 
have checked that they only deform the band structure, 
while preserving the key features summarized above. In 
the numerics on the BHZ model, they are responsible for 
the suppression of the bandwidth of the semi-metal when 
spin-rotational symmetry is broken. 

2. Effective model for the 2D GB inside a 3D TBI 

Similar analysis can be carried out to derive an effec¬ 
tive minimal model for the 2D GB. Let us assume that 
the dislocation sheets are inserted in ^ — 2 ;-plane, such 
that the dislocation lines run in 2 ;-direction, while the 
GB is oriented in x-direction. The coupled surface states 
in the presence of the dislocation sheet are then described 
by the Hamiltonian 

Po = v{kzcr^ + kyCFy)y^ + m{y)y!^, (B4) 

which is a direct generalization of the 2D case, kz is 
still a good quantum number and soliton solutions at 
the dislocation can be found immediately considering the 
kz = ^ case. The projector onto their subspace is now 
given by P = (1 + cr^/i^)/2. 

Allowing again for a generic spin texture along the hy¬ 
bridization direction and including also the inversion and 
and TRS breaking terms, in direct analogy to the 2D 
case, we arrive at the effective model 

= m • V + vkzh^ + (h • v + t sing)d^, (B5) 

where we have defined the effective spin operators = 
a^y^ and = y^. For v = (0,1, 0) this reduces again to 
the case presented in the main text. As with the 2D case, 
we find that the key property determining the response 
to symmetry breaking perturbations is the spin texture 
of the surface states. 


a. Fermi arc-like surface states on GB edges 

In the presence of the inversion symmmetry (m-v=0), 
the effective model possesses also chiral symmetry. This 
means that it can be viewed as a ID SSH model [20] . 
but with a GB momentum g dependent mass M(g) = 
h • V -1-1 sin g. Due to the chiral symmetry, this model has 



two topologically distinct phases: for M{q) < 0 the band 
structure is topologically trivial and characterized by a 
winding number z/ = 0, while for M{q) >0 the winding 
number evaluates to = ±1 and the system is in a topo¬ 
logical phase and exhibits edge states. Thus when the 
TRS breaking is negligible compared to the hybridiza¬ 
tion, the mass depends on the momentum q along the 
GB such that M{q) > 0 for 0 < g' < tt, while M{q) < 0 
for TT < g < 27r, and the edge states of the GB semi-metal 
span only half of the ID edge Brillouin zone of the GB. 

Thus when the 2D GB inside a 3D TBI terminates 
on the 2D surface, we predict the extended GB state to 
contribute additional Fermi arc-like features on the sur¬ 
face states. These states appear as localized flat bands 
at the GB edges parallel to the Burgers vector describ¬ 
ing the dislocations that constitute the GB. Since in a 
real material the GB edge BZ spans only the reduced BZ 
0 < q < 27r/(i, these features are expected to be more 
pronounced for GB defects with a large opening angle. 
As long as they do not hybridize with the TBI surface 
states, which does not happen as long as the surface cones 
occur at different momenta, they can co-exist and could 
in principle be detected via low-frequency optical conduc¬ 
tivity measurements m- Alternatively, breaking inver¬ 
sion symmetry will in general shift the cones to different 
energies, which makes the flat edge states dispersing and 
can lead to edge currents at the GB edges m- 


Appendix C: ID tight-binding picture for spinon 
hy br idizat ion 


Here we present a complementery picture in terms of 
localized spinon wavefunction hybridization. It enables 
us to directly connect the spinon hopping magnitude t to 
the hybridization energy e that arises when two spinon 
binding dislocations are in proximity of each other. 

In the context of the 2D BHZ model, the localized 
spinon zero modes can be solved an¬ 

alytically for a vanishing Rashba coupling near M/B ^ S 
[laiiT]. For single spin component with the dislocation 
located at the origin, the modes are given by the orbital 
space spinors 




sin(Ar) 

\ i 


(Cl) 


where we use polar coordinates r = (r, 6>), ^ = 2B is 
the correlation length, A = Vl — AMB/2B and is a 
normalisation constant. TRS requires that 04,(r) = 0|(r) 
and guarantees that even if these two modes are localised 
on the same dislocation, they do not hybridize even in the 
presence of a Rashba coupling (although the wavefunc- 
tions and the orientation of the spinors in the combined 
spin-orbital space will be modified). Thus isolated dislo¬ 
cations deep in the M-phase bind a Kramers pair of ex¬ 
ponentially localized modes, whose wavefunctions exhibit 
spatial oscillations at the wavelength A, which depends 
on the microscopic parameters. 


In the presence of many dislocations, all the modes 
bound to each dislocation remain near zero energy 
as long as the separation d between all the dislocations 
satisfies d ^ When two dislocations are brought into 
proximity (d the wavefunctions overlap. This gives 
rise to a finite (possibly spin-flipping) tunneling ampli¬ 


tudes ^ f dr^J ^(r)05/^j(r). This hybridizes the 

spinon modes and gives rise to a pair of Kramers degen¬ 
erate modes at finite energy e ^ t (assuming Sz conser¬ 
vation, i.e. t = = til = 0)5 localized 

simultaneously on both dislocations. In principle, the hy¬ 
bridization energy e could be analytically obtained, but 
as the analytic solution is available only in the M/B 8 
regime, we perform instead a numerical analysis. Figure 
shows that the two-dislocation hybridization energy 
obeys the qualitative form e(d) (x cos(Ad)e“^/‘^ through¬ 
out the M-phase, with the wavelength A decreasing with 
smaller M/B. These oscillations in the hybridization en¬ 
ergy, that can cause change the sign as the function of 
dislocation separation d and thus energetically favour dif¬ 
ferent specific separations, follow directly from the oscil¬ 
lating tails of (Gl) and they are qualitatively similar to 
TT-flux vortices in other Dirac-like systems [51 ESI. We 
also find that while the localized spinon modes (Gl) pos¬ 
sess full rotational symmetry near M/B 8, for smaller 
values they become anisotropic as the band gap minimum 
switches from the M-point (full rotation symmetry) to 
the X- and K-points (C4 symmetry). It follows that also 
the hybridization energy becomes anisotropic depending 
on the microscopic parameters. This may have conse¬ 
quences on the nature of the emergent state if the several 
GBs formed arrays of same dimensionality as the parent 
state. 


A grain boundary at small opening angles 0 is well ap¬ 
proximated by a ID array of lattice dislocations. These 
again hybridize pairwise, with the hybridization energy 
e as the function of the GB opening angle 0 shown in 
Figure [^. It is thus natural to ask whether the band¬ 
width of the ID GB can be traced back to the pairwise 
hybridization energy e? To study this, we take t = 2e 
as the simplest ansatz. This is motivated by the fact 
that TRS and translationally-invariant spinful free ID 
fermions with imaginary tunneling amplitude ie have the 
spectrum E± = ±2esing. Figure shows that as the 
opening angle 0 of the grain boundary gets smaller, i.e. 
the spacing d = |b|/ tan^ between the dislocations be¬ 
comes larger and nearest-neighbour approximation is jus¬ 
tified, the spectrum of the mig-gap GB bands in the full 
BHZ model is indeed accurately approximated with the 
relevant two-spinon hybridization energy e{d) as the only 
input. This directly connects the more general edge state 
construction to the hybridization of the localized spinon 
modes and corroborates this picture in the context of the 
BHZ models. 







• Full model 

FIG. 5. A The hybridization energy e as the function of the separation d between two dislocations in units of the bulk band 
gap A. For different values of M/B the hybridization energy always exhibits exponential decay and characteristic oscillations 
that exhibit clear sinusoidal oscillations at larger separations (inset). The kinks as functions of d are artefacts of simulating 
smooth dislocation transport by adiabatically switching on and off the relevant couplings. B The two dislocation hybridization 
energy e in the GB geometry as the function of the GB angle 0 in units of the bulk band gap A. C The mid-gap spectra 
of the 2D BHZ model across the reduced Brillouin zone 0 < kx < 27r/d for different GB angles 0 = arctan(|6|/d) and the 
approximation of each by the unperturbed minimal effective model H 2 D = t sin q over 0 < g < 27r with the hybridization energy 
e{d) = t/2 as the only input. As the dislocation separation d increases, the nearest-neighbour approximation becomes 

more valid and we find excellent agreement. All data has been produced using a 90 x 60 system with M/B — 6. 


Appendix D: Disorder analysis of the grain 
boundary states 

Here we present numerical data on the stability of the 
graphene-like spinon semi-metal when random chemical 
potential disorder is included in the parent TBI. The GB 
states are expected to be stable for two physically moti¬ 
vated complementary reasons: (i) The spinous are rooted 
in the edge states, which can not be removed by moderate 
disorder and (ii) the localized spinon wavefunctions are 
only locally deformed, which can only result in random 


tunneling along the GB, but not in a random chemical 
potential that can gap the modes out. From the point of 
view of the spinon tight-binding model the latter trans¬ 
lates to a random gauge field around the chiral cones that 
only shifts them around, but can not gap them out unless 
the mean displacement is of the order of the cone separa¬ 
tion in momentum space. Thus large opening angle GBs 
are expected to be most stable. 

To substantiate these general arguments, we have stud¬ 
ied the effect of local random disorder on the bi-partite 
hybridization energies e that enter the ID tight-binding 
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FIG. 6. Stability of the graphene-like spinon semi-metal. A Inclusion of a Rashba term of magnitude R breaks Sz conservation, 
which leads to suppression of the bandwidth and different velocities at the cones, but TRS guarantees that the cones are stable. 
B An Sz conservation preserving Zeeman term of magnitude hz, gives rise to a mass for the spinons. In the hybridized spectrum 
it shifts the spinon bands to opposite directions in energy and thereby moves the cones towards each other. The cones annihilate 
and gap out at a quadratic band touching point when hz — ‘le. C When local chemical potential disorder of magnitude w is 
introduced, i.e. /x is a random variable picked from box distribution [(1 — w)(M — 4B), (1 + w){M — 4B)], the hybridization 
energies e become also disordered with the fluctuations, as quantified by the standard deviation Std(e) around the mean (e) 
increasing monotonously. The rate of increase depends on the microscopic details around the nearest-neighbour two dislocation 
hybridization energy corresponding to different GB angles 0 as show in Figure [^. D In the full 2D TBI, we find that the 
gapless spinon semi-metal persists at least up to 50% chemical potential disorder (n enumerates the states), which is insufficient 
to close the bulk gap (inset). We have also simulated randomized dislocation positions (dislocations along a GB of spacing d 
are randomly displaced by up to d lattice constants) and find similar stability. All data has been produced using a 90 x 60 
system in 2D with M/B — 6. The disorder data has been produced for system size 42 x 30 with periodic boundary conditions 
and averaged over 100 disorder realizations. 


model as tunneling amplitudes. Figureshows data for 
the 2D BHZ model demonstrating that even if the dis¬ 
order averaged hybridization energies are strongly influ¬ 
enced by local random disorder, disorder strengths of up 
to 50% are insufficient to cause random sign flips. Thus 
while any disorder potential contains a Fourier compo¬ 
nent that corresponds to nodal scattering (perfect sign 
dimerization would gap out nodes even at maximal sepa¬ 
ration of tt), we expect the spinon semi-metals to be very 
stable. Indeed, numerical data presented in Figure 
shows remarkable resilience with gapless mid-gap bands 
again persisting for disorder strengths of up to 50%. Due 
to the general arguments given above, similar stability 
applies also to the 2D GB semi-metal in the 3D BHZ 


model. 

This stability analysis also implies that perfect trans¬ 
lational invariance of the GB is neither a necessary con¬ 
dition for the stability of the GB semi-metals. Ran¬ 
dom spinon tunneling can arise also due to random turns 
along the GB and thus the stability under local random 
disorder also implies stability in the absence of disor¬ 
der, but with broken translational invariance along the 
GB. Indeed, Figure shows that the gapless mid-gap 
state persists also for randomized GB positions along the 
GB. We therefore conclude that translational invariance 
on average is sufficient for graphene-like semi-metal to 
emerge along them. 

Finally, in Figure we have also considered effects 



















of TRS breaking Zeeman terms. In agreement with the 
effective minimal model derived from the edge states, we 
find hz shifting the Sz conserving spectrum in energy 
and gapping it when hz > 2t. No other Zeeman term 
orientation affects the spectrum as the edge spins are 
polarized in ^-direction. When Rashba term is included 
and the edge states pick up a component parallel to the 


propagation in ^-direction, we find also hy having similar 
effect. Likewise, we have checked the same edge state 
spin texture dependence of the inversion breaking terms 
that, consistent again with the effective model, in general 
shift the nodes to different energies, but do not gap them 
out. 
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